Analog model for quantum gravity effects: phonons in random fluids 
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We describe an analog model for quantum gravity effects in condensed matter physics. The 
situation discussed is that of phonons propagating in a fluid with a random velocity wave equation. 
We consider that there are random fluctuations in the reciprocal of the bulk modulus of the system 
and study free phonons in the presence of Gaussian colored noise with zero mean. We show that in 
this model, after performing the random averages over the noise function a free conventional scalar 
quantum field theory describing free phonons becomes a self-interacting model. 
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Unruh has shown that the propagation of sound waves 
in a hypersonic fluid is equivalent to the propagation of 
scalar waves in a black hole spacetime [l| . Quantizing the 
acoustic wave in such physical system with sonic horizon 
implies that the sonic black hole can emit sound waves 
with a thermal spectrum, recovering the Hawking result 
Q in this analog model, i.e. the presence of phononic 
Hawking radiation from the acoustic horizon. Since the 
Unruh fundamental paper, the possibility of simulating 
aspects of general relativity and quantum fields in curved 
spacetime through analog models has been widely dis- 
cussed in the literature [3, 4). In the context of quantum 
gravity, most of the work on analogue models has fo- 
cussed attention on violation of Lorentz invariance. The 
key point of the discussion is the expectation of a low- 
energy manifestation of a possible spacetime discreteness 
at the Plank scale. Such a manifestation is precisely vi- 
olation of Lorentz invariance. In these analog models 
interesting dispersion relations are derived Q . 

The aim of the present paper is to take a somewhat 
different look at the way an analog model for quantum 
gravity effects can be implemented. Here we have in mind 
the results obtained by Ford and collaborators [7[ , where 
fluctuations of the geometry of spacetime caused by a 
bath of squeezed gravitons has the effect of smearing out 
the light cone. In addition to quantum mechanical metric 
fluctuations, there are induced metric fluctuations gener- 
ated by quantum fluctuations of matter fields. In the 
regime where induced metric fluctuations dominate, the 
conventional approach is to assume a stochastic ensem- 
ble of fluctuating geometries. In this framework, Hu and 
Shiokawa Q assumed that the metric tensor has a deter- 
ministic and a stochastic contribution [9( . We argue that 
in condensed matter physics there is an analog model 
for fluctuations on the light cone generated by a thermal 
bath of squeezed gravitons. We consider a disordered 
medium with random classical fluctuations in the recip- 



rocal of the bulk modulus. Actually, the analog model we 
discuss in the paper is closely related to the induced met- 
ric fluctuations framework. We consider a scalar quan- 
tum field theory associated to acoustic waves defined in a 
fluid with classical random fluctuations, where the veloc- 
ity of propagation of the acoustic waves fluctuates ran- 
domly, and consequently the sonic cone fluctuates. 

Fluctuations of a fluid consist of two types, classi- 
cal and quantum. Quantum fluctuations are due to 
quantization of fluid thermodynamic observables. Con- 
sideration of quantum fluctuations motivated Ford and 
Svaiter [10( to discuss the scattering of light by the zero- 
point density fluctuations of a phonon field. Also, since 
there is an analog to the Casimir effect when the zero- 
point fluctuations of the phonon field are submitted to 
boundary conditions the same authors considered 
in Ref. [12| the effects of boundaries discussing the den- 
sity correlation function associated to the phonon field in 
many different situations. In the present paper we will 
show that fluids are not only useful analog models for 
effects in field theory, but random fluids with particular 
properties can simulate quantum gravity effects. 

Let us start from the basic equations of fluid dynamics 
[l3| . The local mass density p(t, r), the pressure p(t, r), 
the local velocity v(t, r) and the local temperature T(t, r) 
are the thermodynamic fields of the system. For an ideal 
fluid, assuming thermodynamical equilibrium, the acous- 
tic wave equations are obtained by linearizing the fluid 
dynamics equations for small disturbances around the 
fluid at rest, namely 



p(t,r) = po + 5p{t,f), 
p{t,r) = po + Sp(t,r), 
v(t,r) = 5v(t,r), 



(1) 
(2) 
(3) 



where po and po are the constant equilibrium density and 
pressure. Using the Euler and mass balance equations, 
one obtains 
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^Sv(t,r)) + -W5p(t,r)^0, 
dt po 



d_ 

dt 



5p(t,r) + p V ■ Sv(t,f) = 0. 



(4) 
(5) 
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The bulk modulus K of the fluid is defined as K = 
podp/dp — poll 2 , where u = (dp/dpo) 1 ^ 2 is the sound 
velocity. Using that the sound wave in a ideal fluid is adi- 
abatic, 8p(t,r) / 8p{t,r) = dp(t,r)/dpo, one obtains from 
the fluid equations a wave equation for ^{t, r) = Sp(t, r) 

{hw-*)w>* = - (6) 

Here, in general u can depend on r. This wave equation, 
as can be easily checked, is similar as the wave equation 
one would obtain if considering a curved spacetime where 
the metric contains a deterministic and a stochastic con- 
tribution, similar to the situation discussed in Ref. Q, 
where metric fluctuations were represented as fluctua- 
tions in the optical index associated to electromagnetic 
waves. This equation also appears in the study of elemen- 
tary excitations in disordered elastic media with random 
fluctuations, which has been discussed recently by many 
authors. For example, it has been pointed out that a 
free phonon field in a disordered elastic medium becomes 
self-interacting 14]. As we will show in the following, a 
similar situation occurs when u{r) fluctuates. 

When u is a constant, one can derive a wave equation 
also for 8p(t, r) and S(j>(t, r), with Sv = V8<j). In this case 
the quantization of the system follows as usual [l5j ]. the 
phonons have a linear dispersion relation to(k) — u\k\ 
and the operators Sp(t,r) and 5cf>(t,r) obey equal-time 
commutation relation given by [we assume H = 1] 

[Scf>(t,7 f ),8p(t,f')} = -i5 3 (f-?'). (7) 

The unperturbed two-point Green's function of the 
phonon field D F {x,x') is given by the expectation value 
of the time ordered product 

D F (x,x') = -i(0\T(Sp(t,-r)Sp(t',f'))\0), (8) 

where |0) is the ground state of the phonon field. 

Having performed a linearization to obtain a wave 
equation for the acoustic perturbation, let us consider the 
influence of random fluctuations over the phonon field. 
Under such circumstances, the fluid equations become 
random differential equations. The simplest assumption 
one can make is to consider a quantized phonon field in 
a fluid with bulk modulus that fluctuates randomly, with 
Po kept constant (ljj . This generates fluctuations in the 
velocity of the acoustic perturbations. Therefore, let us 
assume (for simplicity we shall consider one-dimensional 
fluids, the extension to three dimensions is not difficult) 

^=4 (1+,/W) ' (9) 

with Mo the sound speed in a homogenous medium and 
v(z) is taken to be a Gaussian colored noise with zero 
mean 

(!/(«)>„ = 0, {v{z) V {z')) v = a 2 h K (z-z'). (10) 



The symbol (... ) v denotes average over noise realiza- 
tions. Under such circumstance, the fluid equations be- 
come random differential equations and two important 
points should be noted. First, one needs to specify a 
stochastic modeling of the fluctuations of the medium 
parameters in terms of random processes. Second, the 
scale separation in the problem should be properly ad- 
dressed. This second point will be discussed latter ahead 
in the paper. With respect to the first point, we mention 
that the specific form of h\ is not relevant for our argu- 
ments, it is sufficient to note that A ~ I /Id, where Id is 
the correlation length of the noise field, and for A — >• oo 
one has uncorrelated white noise. However, it should be 
clear that v(z) is modeling fluctuations of a spacetime 
background on which a quantum field is propagating and 
it should be related to geometry fluctuations on a more 
fundamental level. For example, two-dimensional geome- 
tries are characterized by the scalar curvature R of the 
Riemann tensor and one could speculate that v(z) should 
be related to R. We will not try to make such a connec- 
tion here because this would take us away from the main 
scope of the paper. 

Our aim is to calculate the density-density correla- 
tion function, which is related to an observable as e.g. 
in a light (or neutron) scattering experiment. This is 
most easily calculated solving first Eq. ([6]) and then us- 
ing the relation ip(t, z) = u 2 (z) Sp(t, z). We solve Eq. ([6]) 
by a conventional perturbation expansion in v(z). We 
write ip(t,z) = J2^=o ^n(t,z), where ip n (t, z) is of or- 
der n in v{z). Substituting this in Eq. (|6]) and collecting 
terms of the same order in v{z), one obtains the recursive 
relation for n > 

For n = 0, we have the standard wave equation for 
■00 (t, z). Choosing the solution of the homogeneous equa- 
tion to be zero, the general solution of Eq. <jTTJ) can be 
written as 

r v(z') d 2 

tl) n (t,z)= / dz'dt' G F (t-t';z-z') ' — 2 ^n-ift 7 , z'), 

(12) 

where GF(t — t';z — z') is the Green's function associated 
with the differential operator on the l.h.s. of Eq. (jlll) . 
Introducing a Fourier representation for tp n (t, z) 

Mt,z) = J ^e-^nCw,*), (13) 

and a similar one for Gp(t — t';z — z'), the complete 
solution can be written as 
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^) / n dz - c ^ 

u/ J m =l 



w; z ro _i - z m )v(z m ) ipo{w, z n ), 



(14) 



where it is to be understood that zq = z. 

Given the solution, one can calculate the density- 
density correlation function using Sp(t, z) — ip(t, z)/u 2 (z) 
and quantize as usual. It is important to call attention to 
the validity of the linear dispersion relation in the case of 
randomness. This situation is similar to the relation of 
the light cone fluctuations to Lorentz symmetry. In the 
same way as light cone fluctuations respect Lorentz sym- 
metry on the average, the linear dispersion relation leads 
to k 2 /uj 2 = (1 + v{z)) k 2 /uj 2 , which should be considered 
as valid on the average. 

Before proceeding, let us abstract for the moment from 
the fact that ip represents the pressure field of the fluid 
and consider it as a scalar field per se. Eq. (|6]) represents 
a wave equation for a free scalar field in a fluctuating 
spacetime metric. Suppose one quantizes the field if) and 
calculate the two-point correlation function by taking the 
statistical and quantum averages ((...)„), i.e. average 
over noise realizations and vacuum expectation value of 
products of field operators. The two-point correlation 
function up to second order in v{z) would the be given 
by (in terms of the Fourier transforms of G F and h\) 

k)1>(u, k')) v ) = {2n) 2 6{u - uj')S{k - k')S^, k), 

(15) 

where 

5,/,(a;, k) — iG F {uj, k) + iG F (cu, k)H(ui, k)iG F (u), k), 

(16) 

with 

4 f 

n(«, k) = -2a 2 — / dqiG F {iO, q)h A (k - q), (17) 

where we used the fact that (0\T(ip o (t,z)ip o (t',z'))\0) = 
iG F (t-t',z-z'). 

It is clear that Eq. (flTf leads to the interpretation that 
a non-interacting scalar field propagating in a fluctuat- 
ing metric turns it into a self-interacting field. Notice 



that our results are valid only if the background fluctua- 
tions are much larger than those fluctuations defined in 
equations (HJ-Q. Nevertheless, the scale separation in 
the problem mentioned above can be analyzed in terms 
of its low- and high-frequency regimes. We notice that 
the contributions of the perturbations to the two-point 
(and n-point) functions have a dependence on u) in a 
multiplicative way. More precisely, the deviation from a 
uniform medium occurs multiplicatively with frequency. 
The consequences of this observation are clear; for large 
frequencies, the contributions of the perturbations dom- 
inate. In the opposite regime, for small frequencies, the 
acoustic waves scatter weakly, which implies that, for 
uj — > 0, the field remains essentially non-interacting. At 
this point, one may wonder where the correlation length 
Id fits in this whole scenario. From Eq. ([T7| . we see 
that, for A — > oo the loop integral is ultraviolet diver- 
gent. So, introducing correlations between the classical 
random fluctuations f(z) leads to regularized outcomes. 
Therefore a physically motivated way to obtain a non- 
divergent result is to use colored noise rather than a white 
noise. Another point is that in principle one could think 
of two different ways of performing the joint quantum and 
noise averages. In a static disordered medium one could 
first perform the average over quantum fluctuations for a 
given frozen disorder configuration and then average over 
the disorder configurations 17| . This would correspond 
in random magnetic systems to the case of a quenched 
disorder. In the fluctuating medium considered in the 
present paper, the quantum and disordered averages can- 
not be separated and must be performed simultaneously 
- see e.g. Ref. (l8j | . This case is exactly the case of an- 
nealed systems, where the impurity degrees of freedom 
are in equilibrium with the others degrees of freedom of 
the system. 

Coming back to fluid dynamics, the density-density 
correlation function can be calculated from the relation 
ip(t,z) — u 2 (z)Sp(t, z). From Eqs. ((9]) and (fT4|). one has 



8p{u,z) 



6p (u, z) + v{z)8p {u, z) + XJ(-l)" i^-J J 

/n 
Y[ dz m D F (u,z - z m )v(z m ) 



Y[ dz m D F (uj,z - z m )v(z m ) + 



Sp (uj,z n ), 



(18) 



where 5po(ui,z) — ipo(u>, z)/uq and we used D F for the 8p(t,z) field since its Green's function obeys the same 
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equation as Gf for the ip(t, z) field. Taking the noise and 
quantum averages, the correlation function up to second 
order in v(z) is given in Fourier space by 

((6p(u, k)S P {w', k')) v ) = (2tt) 2 S{u - uj')6(k - k')S p (uj, k), 

(19) 

with 



S p (uj,k) = S^(cj,k) + a 



dq 



iD F (uj,q) 



— 4i^Dp(uj,k)DF(LU,q) 



h A (k - q). (20) 



Note that the induced interactions are different for the 
5p field from those for the i/j field. This is due to the 
wave equation for 5p is not of the form of Eq. ([5]), since 
on replacing tf;(t, z) = u 2 (z)8p{t, z) in this equation there 
are extra spatial derivatives of u(z). Therefore, although 
similar, the model here is different from the ones studied 
in condensed matter physics [T3|. As before, the colored 
noise acts as a regulator of the loop integrals. Of course, 
there is a natural cut-off in the high frequency regime, 
not related to A, in the situation where the wavelength 
of the acoustic phonon is of the order of the interatomic 
separation li, i.e., uj~ x « U. For < l L the field 
theoretical description of the acoustic excitation is no 
longer valid. 

In summary, we have constructed a condensed matter 
analog model for fluctuations of the geometry of space- 
time due to quantum gravity effects. The model leads 
to the interpretation that such classical background fluc- 
tuations induce effective interactions on free fields. We 
expect that such sound cone fluctuations can be tested 
in suspensions like colloids since excitations of acous- 
tic modes in these environments are described by ran- 
dom wave equations as the one we discussed in this pa- 
per [l9j]. A possible course of action is to remember the 
reference LJ], where it is shown that a non- linear dis- 
persion relation generated by a random wave equation 
induces peaks in the density of states. Another direction 
is the analysis of localization of acoustic waves. Such 
issue is related to the problem of how sound cone fluc- 
tuations in our analog model can change the variation 
of flight of pulses. The study of wave localization using 
diagrammatic perturbation method was presented in the 
Refs. [UH. It is known that in one dimension any dis- 
order is strong enough to induce exponential localization 
of eigenmodes. On the other hand, a random fluid with 
a supersonic acoustic flow is also an analog model that 
opens the possibility to discuss the effect of the fluctua- 
tion of the geometry in the Hawking radiation. Specific 
examples studying Anderson localization in the analog 
model of the present paper as well as the implementation 
of randomness in an acoustic black hole will be reported 
elsewhere. 
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